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Zurek has provided a simple picture for the onset of the A- 
transition in 4 He, not currently supported by vortex density 
experiments. However, we argue that the seemingly similar 
argument by Zurek that superfiow in an annulus of 4 He at a 
quench will be measurable is still valid. 



As the early universe cooled it underwent a series of 
phase transitions, whose inhomogeneities have observable 
consequences. To understand how such transitions occur 
it is necessary to go beyond the methods of equilibrium 
thermal field theory that identified the transitions in the 
first instance. 

In practice, we often know remarkably little about the 
dynamics of quantum field theories. A simple question 
to ask is the following: In principle, the field correlation 
length diverges at a continuous transition. In reality, it 
does not. What happens? Using simple causal arguments 
Kibble [[l]J^] made estimates of this early field ordering, 
because of the implications for astrophysics. 

There are great difficulties in converting predictions for 
the early universe into experimental observations. Zurek 
suggested 0j that similar arguments were applicable to 
condensed matter systems for which direct experiments 
could be performed. In particular, for 4 He he argued 
that the measurement of superfiow at a quench provided 
a simple test of these ideas. We present a brief summary 
of his argument. 

Assume that the dynamics of the A He lambda- 
transition can be derived from an explicitly time- 
dependent Landau- Ginzburg free energy of the form 



F(T) = J d 3 x (|W| 2 + «(T)H 2 + i/3|0| 4 ) , 



(1) 



in which a(T) vanishes at the critical temperature T c . 
Explicitly, let us assume the mean- field result a(T) = 
aoe(T c ), where e = (T/T c — 1), remains valid as T/T c 
varies with time t. In particular, we first take a(t) = 
a(T(t)) — —a$t/TQ in the vicinity of T c . Then the fun- 
damental length and time scales £o and to are given from 
Eq.|I] as £q = h 2 /2mao and ro = h/otQ. It follows that 
the equilibrium correlation length £ e q(t) — £, eq (T(t)) and 
the relaxation time r(t) diverge at T c , which we take to 
be when t vanishes, as 
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Although £eq(t) diverges at t — this is not the case 
for the true correlation length £(t), which can only grow 
so far in a finite time. Initially, for t < 0, when we 
are far from the transition, we can assume that the field 
correlation length £ (t) tracks £ eg (t) approximately. How- 
ever, as we get closer to the transition £ e q(t) begins to 
increase arbitrarily fast. As a crude upper bound, the 
true correlation length fails to keep up with £ eg (t) by 
the time —t at which £ eq is growing at the speed of 
sound c(t) = £eq(t)/T(t), which determines the rate at 
which the order-parameter can change. The condition 
d£ eq (t)/dt = c(t) is satisfied at t = —i, where i — ^tqt$, 
with corresponding correlation length 
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After this time it is assumed that the relaxation time is 
so long that £(t) is essentially frozen in at £ until time 
t w +i, when it sets the scale for the onset of the broken 
phase. 

A concrete realisation of how the freezing sets in is pro- 
vided by the time-dependent Landau- Ginzburg (TDLG) 
equation for F of (@) (§, 
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for <j) = [<pi + ifo)/^/^, where r) a is Gaussian noise. We 
can show self-consistently || that, for the relevant time- 
interval — t < t < t the self-interaction term can be ne- 
glected ((3 = 0), whereby a simple calculation finds £ w £ 
in this interval, as predicted. It thus happens that, at the 
onset of the phase transition, the field fluctuations are 
approximately Gaussian. The field phases e l6 ^ r \ where 
4>(r) = \(f)(v)\e l9 ( r \ are then correlated on the same scale 
as the fields. 

Consider a closed path in the bulk superfluid with cir- 
cumference C ^S> £(t). Naively, the number of 'regions' 
through which this path passes in which the phase is 
correlated is Af = 0(C/£(t)). Assuming an independent 
choice of phase in each 'region', the r.m.s phase difference 
along the path is 



A9 C ~ VJ7 = oiy/c/w)). 
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If we now consider a quench in an annular container 
of similar circumference C of superfluid A He and radius 
I -C C, Zurek suggested that the phase locked in is also 
given by Eq|, with £ of Eq.|. Since the phase gradient is 
directly proportional to the superfiow velocity we expect 
a flow after the quench with r.m.s velocity 
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provided I = 0(£). Although in bulk fluid this superflow 
will disperse, if it is constrained to a narrow annulus it 
should persist, and although not large is measurable. 

In addition to this experiment, Zurek also suggested 
that the same correlation length £ should characterise the 
separation of vortices in a quench. In an earlier paper 
H one of us showed that this is too simple. Causality 
arguments are not enough, and whether vortices form on 
this scale is also determined by the thermal activation of 
the Ginzberg regime, in which all 4 He experiments take 
place. Experimentally, this seems to be the case M. 

Our aim in this paper is to see whether thermal fluc- 
tuations interfere with the prediction Eq|| for which ex- 
periments have yet to be performed. Again consider a 
circular path in the bulk fluid (in the 1-2 plane), circum- 
ference C, the boundary of a surface S. For given field 
configurations <f> a (x) the phase change 9c along the path 
can be expressed as the surface integral 



6 c = 2ir j d 2 x p(x), 
Jxes 



(7) 



where the topological density p(x) is given by 

p(x) = 5 2 [0(x)]e jfc a j 1 (x)9 fe 2 (x), i,j = 1,2 (8) 

where ei2 = —£21 = 1, otherwise zero. 

The ensemble average (p(x)) t is taken to be zero 
at all times t, guaranteed by taking (<p a (x))t = = 
(4> a (x)dj4>b{'x))t- That is, we quench from an initial state 
with no rotation. For the Gaussian fluctuations that are 
relevant for the times of interest |||]], a U correlations 
are given in terms of the diagonal equal-time correlation 
function G(r,t), defined by 



b a { X )MO))t=S ab G(r,t) r=|x|. 



(9) 



The correlation length is defined by G(r,t) = 

o( e -'7«t)) j f or large r > Tne TDLG does not lead 

to simple exponential behaviour, but there is no difficulty 
in defining in practice 

The variance in the phase change around C, A9c is 
determined from 

(A8 C ) 2 =^ 2 [ d 2 x f d 2 y(p(x)p(y)) t . (10) 
Jxes J y eS 

The properties of densities for Gaussian fields have 
been studied in detail @||. Define f(r,t) by f(r,t) = 
G(r, t)/G(0, t). On using the conservation of charge 



d 2 x(p(x)p(0)) t =0 
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it is not difficult to show, from the results of that 
A9c satisfies 



(A9 c ) 2 = -[ d 2 x f d 2 yC(\x-y\,t), 

Jx<?S JyeS 

where x and y are in the plane of S, and 

CM)- 1 *' /,2(M) 
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Since G(r, i) is short-ranged C(r, t) is short-ranged also. 
With x outside S, and y inside S, all the contribution 
to (A6>c) 2 comes from the vicinity of the boundary of S, 
rather than the whole area. That is, if we removed all 
fluid except for a strip from the neighbourhood of the 
contour C we would still have the same result. This sup- 
ports the assertion by Zurek that the correlation length 
for phase variation in bulk fluid is also appropriate for an- 
nular flow. The purpose of the annulus (more exactly, a 
circular capillary of circumference C with radius I <C C) 
is to stop this flow dissipating into the bulk fluid. 

More precisely, suppose that C £(t). Then, if we 
take the width 21 of the strip around the contour to be 
larger than the correlation length of C(r, t), Eq.[l^ can be 
written as 



/>oo 
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The linear dependence on C is purely a result of Gaussian 
fluctuations. 

Insofar as we can identify the bulk correlation with the 
annular correlation, instead of Eq.^J, we have 
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The step length £ s (i) is given by 
1 
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There are two important differences between Eq.|l5] and 
Eq. 6|. The first is in the choice of time for which Av of 
Eq.l5| is to be evaluated. In EqJ| the time is the time 
—i of freezing in of the field correlation. Since £(i) does 
not change much in the interval — t < t < i we can as 
well take t = 0. We shall argue below that for Eqjl^ a 
more appropriate time is the spinodal time t sp at which 
the transition has completed itself in the sense that the 
fields have begun to populate the ground states. 

Secondly, a priori there is no reason to identify £, s (t sp ) 
with either £ (or even £(t S p))- In particular, because £ 
in Eq.^| is defined from the large-distance behaviour of 
G(r, t), and thereby on the position of the nearest sin- 
gularity of G(k,t) in the fc-plane, it does not depend on 
the scale at which we observe the fluid. This is not the 
case for £ s {t) which, from Eq.[l6], explores all distance 
scales. Because of the fractal nature of the short wave- 
length fluctuations, £ s (i) will depend on how many are 
included, i.e. the scale at which we look. If we quench 
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in an annular capillary of radius I much smaller than its 
circumference, we are, essentially, coarsegraining to that 
scale. That is, the observed variance in the flux along 
the annulus is nl 2 Av for Av averaged on a scale I. We 
make the approximation that that is the major effect of 
quenching in an annulus. This cannot be wholly true, 
but it is plausible if the annulus is not too narrow for 
boundary effects to be important. 

Provisionally we introduce a coarsegraining by hand, 
modifying G{r,t) by damping short wavelengths 0(1) as 



G(r,t;l) = / d s ke lkx G(k,t)e 



(17) 



We shall denote the value of £ s obtained from Eq. 17 as 



£s(t; I). It permits an expansion in terms of the moments 
of G(k,t)e- k2l \ 



G n (t;l) = / dkk 2n G(k,t)e 
Jo 
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For small r it follows that f' 2 (r, t; I) /{I - f 2 (r, t; I)) 
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Although, for larger, f'(r,t;l) 2 = o(e~ 2r ^^), we find 
that the bulk of the integral Eq.|l6| lies in the forward 
peak, and that a good upper bound for £ s is given by just 
integrating the quadratic term, whence 
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with the equality slightly overestimated. In units of £o 
and Tq we have, in the linear regime la] , 



G n (t;0 
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where I n = J Q dkk 2n e - k2 . The presence of the T(t')/T c 
term is a reminder that the strength of the noise r] is pro- 
portional to temperature. However, for the time scales 
0(t) <C tq of interest to us this ratio remains near to 
unity and we ignore it. For small relative times the inte- 
grand gets a large contribution from the ultraviolet cutoff 
dependent lower endpoint, increasing as n increases. 

If we return to the Landau-Ginzberg equation Eq.|| we 
find that (\4>\ 2 ) t < a /(3 in the interval -t < t < t. 
Although the field has frozen in, the fluctuations have 
amplitudes that are more or less uniform across all wave- 
lengths. As a result, what we see depends totally on 
the scale at which we look. Specifically, from Eq.|2l| 
(™"(0; I) = 0(Z), as shown in the lowest curve of Fig.l. 

If, as suggested by Zurek, we take I — 0(£) we recover 
Eq|j| qualitatively, although a wider bore would give a 
correspondingly smaller flow. However, this is not the 



time at which to look for superflow since, although the 
field correlation length may have frozen in by t = 0, 
the symmetry breaking has not begun. 

Assuming the linearised |B| Eq.|[for small times t > 
we see that, as the unfreezing occurs, long wavelength 
modes with k 2 < t/TQ grow exponentially and soon begin 
to dominate the correlation functions. How long a time 
we have depends on the self-coupling (3 which, through 
Gi, sets the shortest time scale. This is because, at the 
absolute latest, Gi must stop its exponential growth at 
t = t sp , when (|</>| 2 ) tsp , satisfies (|0| 2 )t sp = oto/0. We fur- 
ther suppose that the effect of the backreaction that stops 
the growth initially freezes in any structure. In Fig.l we 
also show £™ m (t; I) for t = M and t = it, increasing as t 



increases 
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For He with quenches of milliseconds the field magni- 
tude has grown to its equilibrium value before the scale- 
dependence has stopped For vortex formation, for 
which the scale is O(£o), the thickness of a vortex, the 
dependence of the density on scale makes the interpre- 
tation of observations problematic. This is not the same 
here. That the incoherent £ s depends on radius I is im- 
material. The end result is that 



m\ C£ s (t sp ;l) 



(22) 



We saw that the expression Eq.|20j for £ s assumed that 
21 is larger than 
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Otherwise the correlations in the bulk fluid from which 
we want to extract annular behaviour are of longer range 
than the annulus thickness. Numerically, we find that 
£e//(0, 1) — 21 very accurately at t = 0, but that 
£ e ff(t,l) > 21 for all t > 0. A crude way to accomo- 
date this is to cut off the integral Eq.[l4[ With a little 
effort, we see that the effect of this is that ^ m {t sp , I) of 
Eq.|20| is replaced by 

c ax (t sP , i) = z? in (t sp , i)[i - (1 - 4* 2 /£ e// (* sp , o 2 ) 3/2 r\ 

(24) 

greater than £™ m (t sp ,l) and thereby reducing the flow 
velocity for narrower annuli. These are the dashed curves 
in Fig.l. The effect is largest for small radii / < £, for 
which the approximation of trying to read the behaviour 
of annular flow from bulk behaviour is most suspect. A 
more realistic approach for such narrow capillaries is to 
treat the system as one-dimensional For this reason 
we have only considered / > £ in Fig.l. We would expect, 
from Eq.^p], that £ s (t sp ; I) has an upper bound that lies 
somewhere between the curves. 

Once I is very large, so that the power in the fluctu- 
ations is distributed strongly across all wavelengths we 
recover our earlier result, that £ s (i sp ; I) = 0(1). In Fig.l 
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this corresponds to the curves becoming parallel as I in- 
creases for fixed t. However, the change is sufficiently 
slow that annuli, significantly wider than £, for which ex- 
periments are more accessible, will give almost the same 
flow as narrower annuli. This would seem to extend the 
original Zurek prediction of Eq.[| to thicker annuli, de- 
spite our expectations for incoherent flow. However, we 
stress again that caution is necessary, since in the approx- 
imation to characterise an annulus by a coarse-grained 
ring without boundaries we have ignored effects in the 
direction perpendicular to the annulus. In particular, 
the circular cross-section of the tube has not been taken 
into account. One consequence of this is that infinite 
(non-selfintersecting) vortices in the bulk fluid have no 
counterpart in an annulus. Removing such strings will 
have an effect on A6c , since the typical fraction of vor- 
tices in infinite vortices is at the level of 70%. However, 
at the spinodal time the fluctuations in A H e are relatively 
enhanced in the long wavelengths, and such an enhance- 
ment is known to reduce the amount of infinite vortices, 
perhaps to something nearer to 20%. The details of this 
effect (being pursued elsewhere) are unclear but, for the 
sake of argument we take the predictions of the curves in 
Fig. 1 as a rough guide in the vicinity of their minima. 

So far we have avoided the question as to which time 
curves we should follow. This is because t sp itself de- 
pends on the scale I of the spatial volume for which the 
field average achieves its ground state value. In practice 
variation is small, with t sp for A He varying from about 
3i to At as I varies from £ <C £ to I — 10£. Since the 
curves for £ s (t sp ;l) lie so close to one another in Fig.l 
once I > 4£ the scale at which the coarse-grained field 
begins to occupy the ground states becomes largely irrel- 
evant. 

— 1/2 

Since Av only depends on £ s it is not sensitive to 
choice of / > 2£ at the relevant t. Given all these ap- 
proximations our final estimate is (in the cm/sec units of 
Zurek @) 
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for radii of 2£ — 4£, tq of the order of milliseconds and 
C of the order of centimetres, v — 1/2 is the mean-field 
critical exponent above. In principle v should be renor- 
malised to v = 2/3, but the difference to Aw is sufficiently 
small that we shall not bother. Given the uncertainties in 
its derivation the result Eq.^ is indistinguishable from 
Zurek's || (with prefactor 0.4), but for the possibility 
of using somewhat larger annuli. The agreement is, ul- 
timately, one of dimensional analysis, but the coefficient 
could not have been anticipated. How experiments can 
be performed, even with the wider annuli that Eqj2^ and 
Fig.l suggest, is another matter. 

We thank Glykeria Karra, with whom some of this 
work was done. This work is the result of a network 
supported by the European Science Foundation . 
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FIG. L_ £? in (t, I) of Eq.20 (solid lines) plotted against 1, in 
units of £, for t=0, 3£ and 4£. of Eq.24 (dashed 

lines) plotted against 1 for t=0, 3£ and 4£. In each case the 
higher lines correspond to higher values of time 
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